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ABSTRACT 

Context. Using the ISOPHOT Serendipity Survey (ISOSS) at 170 pm a sample of galactic star-forming regions exhibiting very cold 
dust temperatures (< 20 K) and high masses (> 100 M Q ) has been established. 

Aims. We characterise the star-forming content of five regions that were selected as potential sites for early stage high-mass star 
formation. 

Methods. We use SCUBA (JCMT) observations in the submillimeter to identify the dense condensations of cold gas and dust. 
Sensitive mid- to far-infared Spitzer observations with IRAC and MIPS allow us to detect associated young stellar objects. From the 
long-wavelength emission we derive dust temperatures and masses for the identified clumps. A sample of associated mid-infrared 
sources is investigated using infrared color-color diagrams and the comparison to a model SED grid to constrain their evolutionary 
stages and derive estimates for additional parameters like the central mass. 

Results. In every region we identify between one and four submillimeter clumps with projected sizes between 0.1 and 0.4 pc. The 
dust temperatures range from 11.6 to 21.3 K and the estimated clump masses are 2 to 166 M Q . Towards the majority of submillimeter 
peaks we find point sources in the near- to mid-infrared. Most are interpreted as low-mass young stellar objects but we also detect very 
red sources. They probably represent very young and deeply embedded protostars that continue to accrete clump material and may 
reach higher masses. Several candidate intermediate-mass proto- or pre-main-sequence stars embedded in the clumps are identified. 
Conclusions. A subset of four clumps may be massive enough (> 100 M G ) to form high-mass stars and accompanying clusters. The 
absence of stellar precursors with current masses in the high-mass regime leave the type of star formation occuring in the clumps 
unsettled. We confirm the presence of large fractions of cold material as derived from large-scale far-infrared measurements which 
dominates the emission of most clumps and suggests that the star-forming process will continue. 

Key words. ISM: individual objects: ISOSS J19357+1950 - ISM: individual objects: ISOSS J19486+2556 - ISM: individual ob- 
jects: ISOSS J20153+3453 - ISM: individual objects: ISOSS J20298+3559 - ISM: individual objects: ISOSS J22478+6357 - Stars: 
formation 



1. Introduction 

The origin of stars with masses of 1OM or more is a 
lively debated issue. Only recently approaches to establish 
an evolutionary sequence based on obser vations and theoret- 
ical considerations have been presented (Beu ther et al.l [2007; 
Zinn ecker & Yorkel I2007I) . However, the observational basis 
concerning the earliest phases in the formation process of high- 
mass stars remains sparse, as many studies target luminous 
infrared sources revealing embedded stellar precursors asso- 
ciated with hot as well as high-density gas. Younger objects 
that are sometimes addressed as infrared-quiet massive cores 
should be best detected using far-infrared and submillimeter 
dust continuum surveys as they are assumed to contain large 
amounts of cold material. Furthermore, the identification of ex- 
tinction features against the galactic mid-infrared background 
provides a number of candidate objects commonly referred to 
as infrared dark clouds (IRDCs). Only a relatively small num- 
ber of can didate sources have been in v estigated in detail (re- 
cently e.g. iBeuther & Steinackerl I2007I: iRathborne etail l2007t 
Birk mann et al.l2007l) but also results of larger-scale studies sug- 
gest that the lifetim es of massive prestellar cores are very short 
(Mott e et alj|2007b . 
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The thermal emission from cold dust (< 20 K) peaks in 
the wavelength range beyond the IRAS 100 / mi limit. The 
170yum ISOPHOT Serendipity Survey (ISO SS) dLemke et al.l 
1996: iBogun et al]|l996t IStickel et al.ll2007t) carried out during 
the ISO mission lKessleretalJ ll996) provides the first large- 
scale survey in this range and is well-suited to search for ex- 
tremely young and massive star-forming regions. Selecting com- 
pact 170yum sources (FWH M < 3.5') that are asso ciated with an 
IRAS lOOyurn point source feeichma n et al.|[l988l) and molec- 
ular gas emission, a sample of more than 50 massive candidate 
sources has been identified (see lKrausd2003l:[Rrause et al.l2 004, 
and references therein). The large-scale average dust tempera- 
tures derived from the flux ratios at 170 and 100 /mi are about 
18 K or lower and the resulting mass-luminosity ratios are M/L ~ 
0.6 M H / L , implying early sta ges in the process of star forma- 
tion (cf. ISridharan et al]|2002l) . In order to explore the physi- 
cal conditions in this unique sample of cold and massive star- 
forming regions and to search for possible high-mass prestellar 
cores, a multi-wavelength follow-up survey has been launched. 
The presence of very young massive clumps and c ores has been 
confir med for the regions I SOSS J20298+3559 dKrause et al.l 
2003), ISOSS J18364-0221 dBirkmann et alJ 120061) and ISOSS 
J23053+5953 dBirkmann et al.ll2007l) . In this paper we report on 
observations of four more ISOSS star-forming regions and ex- 
tend the study of ISOSS J20298+3559. 
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Our continuum observations with high spatial resolution 
covering the submillimeter to near-infrared wavelength range al- 
low us to characterize the detailed star-forming content. In the 
submillimeter we trace the thermal dust emission and identify 
compact condensations. To constrain their evolutionary stage we 
have conducted mid- and far-infrared observations using Spitzer. 
The long-wavelength data in the far-infrared and submillimeter 
is used to estimate the dust temperatures and total masses of the 
clumps. The sensitive mid-infrared data also reveal very young 
associated stellar objects. Additional deep near-infrared images 
provide a good assessment of source confusion. For the detected 
young stellar objects we perform an evolutionary classification 
from their infrared colors. In the Discussion section we further 
characterize a subset using radiative transfer models. 



2. Observations, Data Reduction and Analysis 

2. 1 . Submillimeter continuum emission 

Submillimeter continuum jiggle maps at 450 pm and 850/im 
were obtained with SCUBA ( Holland et alj Il999l) at the 
James Clerk Maxwell Telescope (JCMT) in July 2001 and 
May 2003 under good atmospheric transmission conditions 
(T850um :$ 0-2). Reduction was done u sing the ORAC-DR 
dJenness & Economoul U999) and SURF ( Jenness & Lightfoot 
119981) software. Photometric calibration is based on maps 
of Uranus acquired shortly before or after the observations . 
Further analysis was done using the MIRIAD (ISault et alii 19951) 
software package and fol lowed the procedure described in 
ISandell & Weintrau3<|2001|): To account for the deviations of the 
JCMT beam from a single Gaussian we used the Uranus maps 
to construct symmetric beam models and deconvolve the maps 
of the target regions. The derived beam sizes are 7.8 - 8.8" at 
450/im and 14.7 - 15.1" at 850/im. The maps are then restored 
using a gaussian of 8" and 14", respectively. By fitting gaussian 
components to the restored maps we first derive submillimeter 
continuum fluxes and deconvolved source sizes at 450 pm. We 
then use these sizes convolved with the 14" gaussian to extract 
850/mi fluxes in order to include emission from the same re- 
gions only. The photometric accuracy is estimated to be 30% at 
450 pm and 20% at 850 pm and the pointing uncertainty of the 
submillimeter maps is 2" rms. 



MIPS dRieke et al.l 120041) imaging at 24 pm and 10 pm and 
the MIPS spectral energy distribution (SED) mode. For the 
imaging observations the basic flux calibrated data (BCD) of 
the Spitzer Science Center (SSC) pipeline was used for fur- 
ther data reduction and analysis. The calibration uncerta inties 
of the data are a bout 2% for IRAC dReach et all 120051) . 4% 
for MI PS 24 ami jE ngelbrachl£ialJ|2QQ2b and 10% for MIPS 
10 pm. dGordon et al.l 120071) . Cosmetic corrections and astro- 
metric refinement was per formed with the MOPEX software 
(Mak ovoz & Marleaul2005l) and final images were combined us- 
ing scripts in IRAF. Aperture photometry and PSF fitting was 
done with the aperture corrections given on the SSC websit^- 
The MIPS SEP mod e calibration is based on a spectrum of a 
Boo dLow et al.l2005l) and the measured MIPS 70 pm fluxes. The 
resulting photometric accuracy is estimated to be 5% (IRAC), 
10% (MIPS 24), and 20% (MIPS 70 and SED). 



3. Observational results: Submillimeter emission 
morphology and associated mid-infrared sources 

The five ISOSS regions are listed in Table Q] and displayed in 
Figs. [T]to [5] In the following we describe their individual sub- 
millimeter morphologies and point out the associated 24 pm 
sources. On the IRAC and near-infrared maps those sources for 
which singular counterparts can be identified at shorter wave- 
lengths are marked. 

3.1. ISOSS J1 9357+1 950 

In this region three submillimeter emission components are re- 
solved at 450/im, two adjacent parts (SMM1 North and South) 
and a south-western part (SMM2) (Fig. Q]). They have decon- 
volved FWHM sizes of about 19", 16" and 20" which corre- 
spond to 0.3 1 -0.39 pc. An arc-shaped extended emission feature 
that stretches from South-east to North-west is visible at 24 pm. 
Two 24 pm point sources are detected associated with SMM1 
North, for SMM1 South and SMM2 one and two sources are 
present respectively. Both SMM1 North and South are also coin- 
cident with emission at 70 pm but are not resolved individually. 
At the location of SMM2 no distinct 70 pm feature is detected. 



2.2. Near-infrared emission 

Near-infrared images in J, H and Ks were taken with the Calar 
Alto 3.5 m telescope using the two prime focus wide field 
cameras Omega2000 dBaumeister et aLll2003l) and OmegaPrime 
dBizenberger et al.l il998). Omega2000 features a field of view 
(FOV) of 1 5 .4 x 1 5 .4 arcmin 2 with a pixel scale of 0.4496" pix~ 1 , 
while the FOV for OmegaPrime is 6.8 x 6.8 arcmin 2 with a pixel 
scale of ~ 0.4" pix -1 . The exposure time in the broad band fil- 
ters was 20 minutes each. The exposures were dithered on source 
to allow for sky su btraction. T he reduction and photometry was 
done using IRAF dTodyll 19931) and GAIAQ and the photometric 
calibration for J, H, and Ks is based on the 2MASS point source 
catalog. 

2.3. Mid- and far-infrared emission 

The Spitzer dWerner et al.l 120041) observations include IRAC 
dFazio et al.l |2004l) imaging in all four photometric bands, 



3.2. ISOSS J1 9486+2556 

Three submillimeter clumps are detected in this region: They are 
located along a chain from north-east to south-west and termed 
SMM1 to SMM3 (Fig.©. All three appear compact in the sub- 
millimeter: The deconvolved FWHM diameters are about 16" 
(SMM1), 10" (SMM2) and 15" (SMM3) at 450/im. The cor- 
responding length scales are 0.14 - 0.22 pc. Both SMM2 and 
SMM3 are associated with bright 24 and 70 pm sources. In the 
case of SMM2, several objects are detected at 24pm and the 
brightest one coincides with the submillimeter peak as well as 
with the 10pm source. For SMM3, one source at 24 and 10pm 
is located at the submillimeter peak. In addition, two areas of ex- 
tended 24 and 70 pm emission are found, one in the south-east 
of SMM2 and a second one towards the southern rim of SMM3. 
The first one is not associated with submillimeter emission. Only 
a faint source at 24 pm is detected towards the center of SMM1 
and it is not associated with compact 70 pm emission. 



1 http://star-www.dur.ac.uk/~pdraper/gaia/gaia.html 



2 http://ssc.spitzer.caltech.edu 
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Table 1. Observed ISOSS star-forming regions. 



Region 


Assoc. IRAS 


RA (12000) 


DEC (12000) 


Distance^ 


1 FIR colour 


ISOSS... 


point source 






in kpc 


inK 


119357+1950 


19335+1944 


19:35:45.9 


+ 19:50:58 


4.0 


17.5 


119486+2556 


19465+2549 


19:48:36.8 


+25:56:55 


2.9 


18.0 


120153+3453 


20134+3444 


20:15:20.9 


+34:53:53 


2.0 


19.0 


120298+3559 


20278+3549 


20:29:48.3 


+35:59:24 


1.8 


16.0 


122478+6357 


22460+6341 


22:47:54.1 


+63:57:11 


4.1 


15.4 



" Kinematic distan ces obtained u sing the galactic rotation model of Brand & Blitz ( 1993) and identification of kinematically associated molec- 
ular cloud complexes (Rrause 2003). 
* T F i Rco i or is derived from the far-infrared flux ratios at 100 fim (IRAS) and 170 /urn (ISOSS). 




48 46 44 42 19:35:40 48 46 44 42 19:35:40 48 46 44 42 19:35:40 

Ri^ht astension 

Fig. 1. Observations of the region ISOSS J19357+1950. The left panel shows the 24/im map and the overlaid contours show the 
SCUBA 450/mi emission. The mid panel shows a color composite of the IRAC images. The right panel is a JHKs color composite 
using OmegaPrime and Omega2000 observations. The circles surround the sources investigated in the Analysis section. 




Right ascension 

Fig. 2. Observations of the region ISOSS J19486+2556. Same arrangement and symbols as in Fig.Q] 



3.3. ISOSS J201 53+3453 

This region contains a single submillimeter clump clearly de- 
tected in both SCUBA bands (Fig. [3]). The deconvolved 450yum 
FWHM extension is about 18" x 14" i.e. the projected diameter 
of the clump is around 0. 16 pc. In the 24 fim band two point-like 
sources surrounded by extended emission are detected towards 
the north-western limb of the submillimeter clump. The brighter 
one is located among a cluster of sources seen at shorter wave- 
lengths. At 70/L/m none of these are resolved individually, they 



blend into a single slightly extended source. The submillimeter 
emission peak is offset from the brighter point source by 11". 

3.4. ISOSS J20298+3559 

This region has been studied in detail in iKrause et al. (2003). 
Four submillimeter emission peaks are found in this region: Two 
of them are joined by extended emission (SMM1 and SMM2) 
while SMM3 and SMM4 are offset to the wes t (Fig. |4j. SMM4 
has not been discussed in IKrause et alj d2003l) . The sizes of the 
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clumps SMM1 and SMM3 are approximately 0.14 and 0.17pc, 
SMM2 is unresolved (< 0.07 pc) and SMM4 is more extended: 
The deconvolved FWHM diameter is roughly 39" correspond- 
ing to 0.34pc. SMM1, SMM2 and SMM4 are associated with 
emission at 70 fim. SMM3 was not covered by our 70 /mi maps. 
Two 24 fim point sources are detected towards SMM2 and one 
towards SMM4. SMM1 and SMM3 have faint counterparts at 
24 fim. 

3.5. ISOSS J22478+6357 

The submillimeter emission in this region traces an elongated 
clump (SMM1) that is resolved in an eastern and a western com- 
ponent at 450/L/m (Pig. [5J. The deconvolved 450/mi FWHM ex- 
tensions of these are ~1" and ~ 12", corresponding to roughly 
0.14 and 0.24 pc. At 24 /mi, several point sources are associ- 
ated. The brightest are located at the eastern and to the southern 
limb of the western submillimeter peaks. The 70 fim map reveals 
emission that coincides with the submillimeter peaks and with 
the 24 fim sources. 



SMM1 




100 1000 
Wavelength in |j.m 



Fig. 8. Spectral energy distribution of the submillimeter source 
in the region ISOSS J20153+3453. The diamonds represent 
photometric data obtained with SCUBA (450 and 850/mi). 
Histogram-like bars show the error range of the MIPS SED spec- 
trophotometry (53-93 fim). The lines show modified blackbody 
fits for A > 93 fim using different dust opacities (solid: no ice 
mantles, dotted: thick ice mantles, dashed: thin ice mantles). 



4. Analysis 

We first present the analysis of the long-wavelength observations 
towards the star-forming regions and then address the associated 
infrared sources. The extracted fluxes in the far-infrared and sub- 
millimeter are listed in Table |2] In Figs. [6] [7] [8] [9] and [10] the 
spectral energy distributions (SEDs) of the detected clumps are 
sketched. 



4.1. Far-infrared spectral slopes 

The examination of the SEDs compiled for the seven detected 
clumps towards which MIPS SED observations have been per- 
formed shows that the spectral slopes at wavelengths around 
70 fim and below do not resemble a single thermal emission 
component that would reproduce the submillimeter fluxes. It 
presumably originates from warmer du st components or tran- 
siently heated very small grains (cf. iBirkmann et alJ l2007t 
iRawlings et al.l [2005b . The spectral slopes (dlog(/iF / |)/dlog/l) 
vary in the range 0.3 - 2.0 for 70 fim < A < 94 fim. 



4.2. Clump dust temperatures and masses 

To characterize the cold component of dust (and gas) that gives 
rise to the emission at long wavelengths, we assume that it can 
be reproduced by isothermal dust emission. As noted above, the 
fluxes at wavelengths around 70 fim and below do not conform 
with this assumption and we therefore use the datapoints at 93, 
450 and 850 fim to estimate the dust temperature. In that respect 
we also assume that the emission at these wavelengths is opti- 
cally thin and can be approximated by a modified Planck spec- 
trum. We use the dust opacities given in lOssenkopf & Hen ning 
(1 1994 for a MRN distribution either without ice mantles, with 
thin and with thick ice mantles (/fsso/jm = 0.6 - 1 .4 cm 2 /g) corre- 
sponding to emissivity indices of f3 — 1 .8 1 - 1 .85 at long wave- 
lengths. The dust-to-gas mass ratio is assumed to be 1/100. In the 
cases where no MIPS SED measurements are available, we ex- 
trapolate the 70 fim flux to 93 fim allowing for spectral slopes of 
0.3 and 2.0 (adopting the measured extremes). We list the result- 
ing range of dust temperatures and gas masses for each clump in 
Table |2] The fitted curves are shown in Fig. [6l 171 [8l |9land[T0l 
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Fig. 4. Observations of the region ISOSS J20298+3559. The left panel shows the 24 //m map and the overlaid contours show the 
SCUBA 450/mi emission (for SMM3 the 850yimi emission is shown). The mid panel shows a color composite of the IRAC images. 
The right panel is a JHKs color composite using OmegaPrime and Omega2000 observations. The circles surround the sources 
investigated in the Analysis section. 
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Fig.5. Observations of the region ISOSS J22478+6357. Same arrangement and symbols as in Fig.Q] 



4.3. Associated infrared sources 

We detect a number of sources at 24 jjm and shorter wavelengths 
that are associated with the submillimeter clumps. To evaluate 
their nature and relationship with the submillimeter condensa- 
tions we select the sources marked in the mid and right pan- 
els of Fig. [TJ to [5] for which we can identify singular counter- 
parts. At distances of 2 - 4kpc the resolution achieved in the 
Ks band corresponds to 1600 - 3200 AU and we can only probe 
the confusion of sources on scales larger than these. For sources 
that show no near-infrared counterpart the IRAC resolution at 
3.6 /im which corresponds to 3000 - 6000 AU limits the search 
for multiplicity. In addition we check for optical counterparts 
in Second Digitized Sky Survey DSS2-red maps. Sources with- 
out optical counterpart are best candidates to represent embed- 
ded objects and therefore the association with the molecular and 
dusty clumps is presumably no projection effect. However, the 
actual positions of the sources with respect to the clumps and 
the foreground extinctions are unknown. 

iRobitaille et alj d2006l) present a grid of young stellar ob- 
jects (YSOs) model SEDs which has been used to analyze 
mid-i nfrared sources detected with Spitzer in star- forming re- 
gions dlndebetouw et alj2007tlShepherd et alj2007l) . Utilizing a 



Monte-Carlo radiative transfer code, more than 20 000 2D mod- 
els were computed for ten inclinations each. The different evolu- 
tionary phases in the modell ed parameter space are ordered ac- 
cording to the Class system (jAndre 1994) into stages: Stage 0/1 
(hereafter: Stage I) are very young objects that are deeply em- 
bedded in an extended envelope with a cavity surrounding an ac- 
cretion disk, Stage II contains disk-dominated models and Stage 
III represents star-dominated systems. Via a web interface mea- 
sured fluxes of individual sources can be put int o a model SEP 
fitting routine to constrain source parameters ( Robitaille et al.l 
2007) and we make use of this below in the Discussion section. 
The color-color diag ram [3.6]-[5.8] versus [8]-[24] (Fig. 23 in 
Robitaille et al. (2006)) offers a good distinction between the 
stages . These colors have also been compiled in iReach et"ail 
d2004i) for a number of young stellar objects in the Elephant 
Trunk Nebula wh ich have been clas sified before according to the 
Class system (see lReacTTe t al. 2004, and references therein). For 
a verification of the color assessment we plot th eir data and the 
region s corresponding to the Stages proposed by Robitaille et al. 
( 2006) in the upper left panel of Fig. QT| The foreground extinc- 
tion appears to be small for these sources. The Class 1/0 and I 
sources fall into the corresponding Stage I region with one ex- 
ception. The colors of Class I/II and II sources lie in the Stage 
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Fig. 6. Spectral energy distributions of the submillimeter sources in the region ISOSS J19357+2556. The diamonds represent 
photometric data obtained with MIPS (70 /mi) and SCUBA (450 and 850 /mi) and extrapolated data at 93 /mi for SMM1 South and 
SMM2. Histogram-like bars show the error range of the MIPS SED spectrophotometry (53-93 /mi) for SMM1 North. The lines 
show modified blackbody fits for A > 93 /mi using different dust opacities (solid: no ice mantles, dotted: thick ice mantles, dashed: 
thin ice mantles). For SMM1 South and SMM2 two sets of curves are shown for spectral slopes of 0.3 and 2.0 in 70 < A < 93 fim. 




Wavelength in urn 

Fig. 7. Spectral energy distributions of the submillimeter sources in the region ISOSS J 19486+2556. The diamonds represent photo- 
metric data obtained with SCUBA (450 and 850 fim). Histogram-like bars show the error range of the MIPS SED spectrophotometry 
(53-93 fim), for SMM1 upper limits are plotted with triangles. The lines show modified blackbody fits for A > 93 fim using different 
dust opacities (solid: no ice mantles, dotted: thick ice mantles, dashed: thin ice mantles). 



II region with one exception. The colors of the debris disk ob- 
ject p Pictoris place it on the border between the Stage II and III 
regions. These accordances support the utility of the colors for 
classification and we apply this method to our source samples to 
get an estimate of their evolutionary status. However, due to the 
presence of a significant amount of dust in the regions, an en- 
hanced extinction is expected for our sources. In the remaining 
five panels of Fig. QT|the infrared colors of the sources marked 
(and numbered) in Fig. Q] to [5] are plotted. Sources without opti- 
cal counterpart in DSS2-red maps are marked with crosses. Most 
of the evaluated sources exhibit colors lying in the Stage II re- 
gion or in the vicinity of the Stage I/II transition. Interestingly, a 
few sources show very red colors that are consistent with Stage 
I even if a foreground extinction of Ay > 100 is assumed. 
Exhibiting [8]-[24] > 5, these sources are distinct from the sam- 
ple in the upper left panel. Two red sources, Source 1 in ISOSS 
J19486+2556 and Source 6 in ISOSS J22478+6357, apparently 
show enhanced 5.8 //m emission but still lie within the Stage I 
region. Both are very faint in the near-infrared. 

Several sources are associated with "fuzzy" green features in 
the IRAC composite images and/or extended emission features 
in the near-infrared: Source 2 in ISOSS J19486+2556, Source 2 
in ISOSS J20153+3453, Source 8 in ISOSS J20298+3559 and 
also towards the SMM3 peak in ISOSS J19486+2556, where we 
cannot identify the appropiate counterparts to the detected 24 /mi 



source, and towards the SMM1 peak in ISOSS J20298+3559, 
where the 24 /mi source shows no counterparts at shorter wave- 
lengths. These features likely trace outflow activity from young 
accreting objects which gives rise to shock-excitation of molec- 
ular hydrogen and CO. However, also atomic hydrogen Br a line 
emission can contribute to the IRAC 4.5 /mi flux. 

The most prominent features in the IRAC composite image 
of ISOSS J19486+2556 are two extended red emission features 
east of SMM2 and south of SMM3 which are also visible at the 
MIPS wavelenghts. In the near-infrared very red objects are de- 
tected that are located towards the centers of the features which 
themselves are not traced. Extended low-density dust heated to 
higher temperatures could give rise to the mid- and far-infrared 
emission. 

In the panel of Fig. QT| corresponding to ISOSS 
J20298+3559 we also plot the approximate colors of the 
brightest source in the region, locat ed western o f SMM 1. It 
has been investigated in detail in iRrause et al. I (120031) and 
classified as Herbig B2e star that is accreting from a disk 
which is optically thick. The 8 and 24 //m images seem to show 
non-linear effects due to the brightness of the source which may 
affect the derived colors corresponding to Stage II. 
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Wavelength in urn 

Fig. 9. Spectral energy distributions of the submillimeter sources in the region ISOSS J20298+3559. The diamonds represent 
photometric data obtained with MIPS (70 /mi) and SCUBA (450 and 850 /mi), extrapolated data at 93 /mi for SMM1 and SMM2 
and the datapoint at 1 .2 mm for SMM1 . Histogram-like bars show the error range of the MIPS SED spectrophotometry (53-93 //m) 
for SMM4. The lines show modified blackbody fits for A > 93 /mi using different dust opacities (solid: no ice mantles, dotted: thick 
ice mantles, dashed: thin ice mantles). For SMM1 and SMM2 two sets of curves are shown for spectral slopes of 0.3 and 2.0 in 
70 < A < 93 /on. 

Table 2. Long-wavelength emission and derived properties of the cold component of the detected clumps. 



Region 
ISOSS. . . 


Clump 


70 /an 


Total flux (Jy) 
93 /an 450 /an 


850/mi 


Size (pc) 


Dust temperature (K) 


Gas mass (M G ) 


J19357+1950 


SMM1 N 


3.3 


6.0 


2.1 


0.40 


0.37 


16.6 - 18.7 


54-92 




SMM1 S 


0.57 


0.80- 1.3" 


1.0 


0.21 


0.31 


14.1 - 16.5 


34-66 




SMM2 


0.21 


0.30 - 0.47 a 


1.9 


0.35 


0.39 


12.0- 13.7 


85 - 166 


J19486+2556 


SMM1 


<0.15 


<0.26 


2.2 


0.51 


0.22 


11.6- 12.5 b 


69 - 123 b 




SMM2 


2.7 


4.6 


2.1 


0.38 


0.14 


16.1 - 18.1 


28-49 




SMM3 


7.4 


11 


4.1 


0.59 


0.21 


16.8 - 19.0 


46-79 


J20153+3453 


SMM1 


13 


17 


12 


2.1 


0.16 


15.3 - 17.0 


87 - 149 


J20298+3559 


SMM1 


1.7 


2.5 -3.7 a 


1.3 C 


0.23 c 


0.1 4 C 


15.8- 18.8 d 


6- 13 d 




SMM2 


1.1 


1.6-2.4 a 


0.31 c 


0.091 c 


< 0.07 c 


17.3-21.3 


2-3 




SMM4 


0.41 


0.81 


3.9 


1.0 


0.34 


12.1 - 13.2 


46-80 


J22478+6357 


SMM1 E 


0.80 


1.6 


1.0 


0.35 


0.14 


14.7 - 16.4 


47 - 81 




SMM1 W 


0.72 


1.0-1.6" 


1.8 


0.47 


0.24 


13.4- 15.5 


79- 153 



a Extrapolated values. 

b Derived using upper limit fluxes. 

c Flux values and sizes from dKrause et al]|2003t) . 

d Derived from listed fluxes and 0.075 Jy at 1.2 mm dKrause et all2003h . 



5. Discussion 

5.1. The clump population 

Our simple approach to reproduce the long-wavelength emission 
towards the sample of submillimeter clumps by one component 
neglects any effects of the density and temperature distributions 
within the clumps and the assumption of optically thin emission 
may be violated in the far-infrared if the (column) density is high 
enough. Assuming a density profile that is fiat in the innermost 
region of 100 AU radius, follows a power-law (p oc r~" with a = 
1 .5) in between and is cut-off exponentially at 20 000 AU, an op- 
tical depth of unity is reached at about 100 /mi for a 100 M Q con- 
densation (dust opacities from lDraine & Leel(ll984l) : Steinacker, 
private communication). Therefore we may underestimate the 
total emission at the shortest wavelengths used (93 //m) and also 
the mean dust temperatures. Setting the dust temperature to a 
canonical value of 20 K results in a decrease of the mass by up 
to a factor of about two for the coldest temperatures, but would 
require high optical depths at wavelengths around 100 /mi to be 
consistent with the measured fluxes. However, the required mass 
scales similarly to the optical depth. It seems therefore unlikely 



that this effect causes errors in the mass estimation of more than 
a factor of two. 

The distances to the regions are mainly derived from radial 
veloci ties and carry therefore some uncertainty. In lKrause et al.l 
(120031) the distance to ISOSS J20298+3559 has been constrained 
by an extinction study to 1.8 + 0.3 kpc. We estimate the other 
distances to be reliable within ~ 1 kpc which results in possi- 
ble mass changes by factors of 1.3 to 1.8. Another uncertainty 
factor of the order of two for the mass estimate arises from 
the unknown dust compo s ition, e.g. the dust opacities given 
by [Ossenkopf & Henningl ([1994) are about five times higher 
than those used for the calculation of the model SED grid 
(IWhitnev et al.ll2003l) . It is obvious from the SED diagrams that 
the used dust models in combination with the one-component 
approach does not fit the data well for all sources: In the cases 
of SMM2 in J20298+3559 and SMM1 East in J22478+6357 it 
seems that dust models corresponding to a lower value of the 
emissivity index /? at long wavelengths or the introduction of 
several emission components would improve the fit. However, 
the latter approach would increase the number of variable pa- 
rameters for the fit above the number of datapoints used. From 
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Fig. 11. IRAC and M IPS infrared color-color diagram. The regions where most models are Stage I, II and III according t o 
Robitaille et al.l d2006l) are shaded from dark to light grey. The vectors show an extinction of Ay— 20 (jlndebetouw et alj |2005). 
In the upper left panel young stellar objects corresponding to different evolutionary stages from iReach et al.l (120041) are plotted. 
In the five remaining panels the mid-infrared sources associated with submillimeter emission in the ISOSS star-forming regions 
are drawn (diamonds, numbering as in Figs. [T]-[5]l. The sources without optical counterpart are marked by crosses. In the panel 
corresponding to ISOSS J20298+3559 the approximate colors of the Herbig B2e star are represented by the triangle. 



the available measurements in the submillimeter range we can- 
not place constraints on the dust properties. 

The derived dust temperatures (see Tabled lie between 11.6 
and 21.3 K, approximately in the same range as the far-infrared 
color temperatures from ISO and IRAS fluxes of the whole re- 
gions (see TableQ]). In comparison to other studies, the values are 
consistent with those of cold dust components derived and as- 
sumed for young dense clumps (e. g. Beuthe r & Steinacker|2007t 
ISridharan et"ai1l2005t lEnoch et alj|2007l) ~ 

The estimated masses that reside in these cold dust and gas 
components range from 2 M up to 166 M . Taking the mass es- 
timates at face value, four of the clumps may be massive enough 
(> 100 M ) to form a high-mass star and an accompanying clus- 
ter, if one assumes t hat a fraction of 10% of th e total mass ends 
up in massive stars dZinnecker & Yorkdl2007l) . The star forma- 
tion efficiency is observationally not well constrained on the 
scales of individual clumps. 

From the source sizes in Table [2] volume-averaged den- 
sities have been derived for our clump sample under the 
assumption of spherical symmetry. The resulting densi- 
ties are 4- 10 4 crrT 3 < nn 2 < 10 6 crrT 3 with a median value 
of 2-10 5 cm~ 3 . Comp ared to studies of d ense condensa- 
tions found in IRDCs ( Rathborn e et al.l 120061, m edian density 
3 • 10 4 cm 3 ) and in Cvgnus X dMotte et alj|2007l median den- 



sity 2-10 cm ) the evaluated clumps seem to be more simi- 
lar to the latter sample based on a millimeter continuum survey. 
Both the mentioned samples seem to contain candidate objects 
representing very early stages of massive star formation (e.g. 



massiv e infrared-quiet cores). For the sample of B euther et al.l 
(2002) based on IRAS measurements and containing luminous 
infrared s ources a low e r mea n density of 9 • 10 3 crrT 3 is recal- 
culated in Mo tte et all (|2007). The direct comparison of these 
densities may be affected by systematically different derivation 
methods, but it suggests that the source selection based on longer 
wavelengths gives access to presumably less evolved sources. In 
the present case the 170/vm band measurements constrain the 
average dust temperature to lower values than those that can be 
inferred from the IRAS bands. 



5.2. The nature of the 24 iim sources 

The majority of the clumps appear in a stage where 24 //m 
sources are present. In the following we discuss their ability 
to generate high-mass star progenitors and evaluate selected 
subsamples in more detail. During their formation high-mass 
stars develop high l uminosities: According to the approach of 
iMcKee & Tar] (2003) values above 1OOL are reached already 
when one solar mass has been accreted. The corresponding 
timescale is on the order of several 10 4 years. Therefore young 
high-mass stars should stand out as luminous sources. However, 
due to the reprocessing of the radiation at short wavelengths by 
the surrounding dusty envelopes only a fraction of the luminos- 
ity is released at the infrared wavelengths we probe (here we 
consider A < 24yum). If the source is deeply embedded, we will 
expect its infrared colors to be very red and the emission at op- 
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Fig. 10. Spectral energy distribution of the submillimeter 
sources in the region ISOSS J22478+6357. The diamonds repre- 
sent photometric data obtained with MIPS (70 /an) and SCUBA 
(450 and 850 /mi) and extrapolated data at 93jum for SMM1 
West. Histogram-like bars show the error range of the MIPS 
SED spectrophotometry (53-93 /urn) for SMM1 East. The lines 
show modified blackbody fits for A > 93 /mi using different 
dust opacities (solid: no ice mantles, dotted: thick ice mantles, 
dashed: thin ice mantles). For SMM1 West two sets of curves 
are shown for spectral slopes of 0.3 and 2.0 in 70 < A < 93 /mi. 



tical wavelengths to be not detectable. By combining the classi- 
fication from the colors with an estimate of the total luminosity 
including the far-infrared and submillimeter data we therefore 
get a handle on the probability of a source to be a high-mass 
star precursor: Sources of Stage II or later should generate high 
luminosities to be considered candidate massive young stellar 
objects. 

We consider the sources' locations in the infrared color- 
color space as qualitatively describing their evolutionary stages. 
From the sample of 24 //m sources that has been evaluated in 
the Analysis section the majority is estimated to be compara- 
ble to Stage II models. As the majority of the associated clumps 
does not substantially exceed a total luminosity of 200 L , we 
consider most of these objects as being rather low-mass YSOs. 
Their association with the submillimeter clumps suggests that 
they have formed inside and that the clumps contain density sub- 
structures. 

5.2.1. Very red protostar candidates 

As described above, sources that display very red infrared col- 
ors and sources that are not detected in the shorter-wavelength 
bands are best candidates to represent very young objects that 
may evolve to massive YSOs. Using Fig. [TT] we select Source 
5 in ISOSS J19357+1950, Source 1 in ISOSS J19486+2556 
and Source 2 in ISOSS J20153+3453. The remaining sources 
with [8] -[24] > 5 are not associated with massive submillimeter 



clumps and probably will not be able to accrete enough mate- 
rial to end up as high-mass stars. Furthermore, in iRrause et alj 
d2003l) the submillimeter source SMM1 in ISOSS J20298+3559 
is found to be a candidate Class object. We detect a 24 /mi 
source towards this peak which is not seen at shorter wave- 
lengths. Associated outflow-tracing emission supports the clas- 
sification as young protostellar object. 

Source 5 in ISOSS J19357+1950 is located towards the 
massive clump SMM2. We apply the online model fitter us- 
ing the fluxes including upper limits in the H band, at 70, 450 
and 850/mi. Source 2 found towards the massive clump ISOSS 
J20153+3453 SMM1 exhibits very red colors and is associ- 
ated with a possible outflow. The online model fitter is fed with 
the fluxes including upper limits in the J band, at 70, 450 and 
850/mi. For SMM1 in ISOSS J20298+3559 we input the fluxes 
measured at 24, 70, 450 and 850//m and upper limits in the 
IRAC bands and in the Ks band. The range of possible fore- 
ground extinction is selected to Ay < 500. The fitting results 
for these three sources are listed in Table [3] From the inspec- 
tion of the SED plots we establish a cut-off value of x\ (x 2 P er 
datapoint) to select the models that fit the datapoints within rea- 
sonable margins. We also give the total extinction to the cen- 
tral source by combining the fitted foreground extinction with 
the extinction along the line of sight from the model outer radii 
inwards. These parameters are degenerate to some extent and 
allow models with different envelope sizes to be fitted via com- 
pensating foreground extinctions. The results are consistent with 
the interpretation of the sources as young embedded stellar pre- 
cursors. The observed outflow activity signatures accord with 
high accretion rates for the latter two sources. In the cases of 
ISOSS J19357+1950 SMM2 and ISOSS J20153+3453 SMM1 
the estimated large total masses of around 100M o residing in 
cold components of dust and gas would in principle allow on- 
going accretion and the build-up of high-mass stars. However, 
due to the multiplicity of near-infrared sources that are de- 
tected in the vicinity of ISOSS 20153+3453 SMM1 the further 
evolution of this object is not straightforward. Towards ISOSS 
J19357+1950 SMM2 there also is a second more evolved source 
nearby. In the third case the mass estimate of about 10M Q for 
ISOSS J20298+3559 SMM1 is consistent with the source being 
an intermediate-mass star precursor. 

The remaining very red candidate Stage I source located to- 
wards the clump SMM1 in the region ISOSS J19486+2556 is 
particularly interesting because of the low infrared luminosity 
in comparison to the submillimeter, making it a candidate for a 
young and deeply embedded object. We input the photometric 
data into the online model fitter including upper limits in the Ks 
band, at 70, 450 and 850/mi. We choose a cut-off at xl < 49.6 
to select the acceptable fits. The results for this source are am- 
biguous: Among the selected models are several belonging to 
the Stage I category with central source masses of about 0. 1 M Q 
but also a majority of Stage II models with central source masses 
between 3 and 4.5 M , both with and without envelopes. There 
is no significant difference in the fitted total extinction towards 
the central sources. We therefore cannot assess and constrain the 
parameters of this source. 

5.2.2. Evolved young stellar objects 

We also detect sources whose SEDs are dominated by emission 
at mid-infrared wavelengths which exhibit higher total luminosi- 
ties and may represent more massive YSOs. By comparing the 
spectral energies at 24 /mi and 70/mi we find two candidate 
sources, located towards SMM4 in ISOSS J20298+3559 (Source 
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Table 3. Results of the SED fitting for three very red young stellar objects. 



Region ISOSS. . . 


J19357+1950 


J20153+3453 


J20298+3559 


Source 


5 


2 




Associated clump 


SMM2 


SMM1 


SMM1 


Selected cut-off 


Xt < 14.2 


Xi < 9.3 


Xi < 12.5 


Fitted models 


Stage I 


Stage I 


Stage I 


Total extinction 


A v ~ 85 


A v ~ 100 


A v > 200 


Central mass (M G ) 


0.2-2.1 


0.7-2.3 


0.2-4.6 


Envelope accretion 








rate (M Q /yr) 


1.6 ■ 10~ 6 - 8.9 • 10" 5 


1.5 • 10~ 6 - 3.4 ■ 10~ 4 


4 • 10~ 6 - 7 • 10~ 4 



Table 4. Results of the SED fitting for two evolved young stellar 
objects. 



Region ISOSS. . . 


J20298+3559 


J22478+6357 


Source 


8 


3 


Associated clump 


SMM4 


SMM1 East 


Selected cut-off 


Xl < 3-3 


*d<2-l 


Fitted models 


Stage IP 


Stage II b 


Total extinction 


6 < A v < 9.4 


29 < A v < 41 


Total luminosity (L Q ) 


~ 150 


~ 2000 


Central mass (M Q ) 


3.2-4.9 


6-8.5 


Preferred inclination 


low 




Disk mass (M Q ) 


4 • 10~ 6 - 10-' 


10~ 5 - 10 1 


Accretion rate (M G /yr) 


io-" - 10~ 6 


5 • IO" 11 - 2- 10 s 


System age (yr) 


10 6 - 10 7 


10 6 - 6 • 10 6 



a M disk /M, > 7.7 • 10- 7 
b M dlsk /M, > 9.5 • IO" 7 



8) and SMM1 East i n ISOSS J22478+63 57 (Source 3). We also 
use the fluxes from lKrause etal.1 (120031 IRS 6) for the I, J, H 
and Ks bands as well as upper limits at 70, 450 and 850 fim for 
the former and upper limits in the J band, at 70, 450 and 850 fim 
for the latter source. The fitting results are compiled in Table |4] 
We consider the first source to be an embedded proto- or pre- 
main-sequence star of intermediate mass as also suggested by 
iKrause et all (|2003). It probably has been formed in the SMM4 
clump. The results for the second source are consistent with 
the source representing an evolved YSO of at least intermedi- 
ate mass that is still embedded and presumably has been formed 
within the clump SMM1 East. 

6. Conclusions 

We have analysed multi-wavelength observations of five star- 
forming regions that were identified using the ISOPHOT 
Serendipity Survey at 170^/m. From the discussed results we in- 
fer: 

1 . We found one to four compact (~ 0.2 pc) submillimeter con- 
densations in every region that represent molecular clumps 
containing a cold component of gas and dust. The dust tem- 
perature estimates vary between 1 1.6 and 21.3 K and accord 
with the large-scale color temperatures measured in the far- 
infrared. 

2. The resulting estimated clump masses that reside in these 
cold components range from 2 to 166 M . Four out of twelve 
clumps may be massive enough (> 100 M ) to be promising 
candidate birthplaces for high-mass stars. 

3. We identify multiple associated mid-infrared sources for the 
majority of the clumps suggesting that they embody signifi- 
cant density substructures on smaller scales than probed by 



the submillimeter observations 0.1 pc). Since the emis- 
sion of most clumps is dominated by the cold material we 
expect thriving star formation from further collapse of over- 
densities. 

4. Most of the associated sources are considered as low-mass 
young stellar objects in an evolutionary state later than Class 
0/1 partly embedded in the clumps where they have formed. 

5. In the region ISOSS J19357+1950 we find a massive clump 
of around 100 M . A probably deeply embedded accret- 
ing Class 0/1 protostar of ~ 1 M is detected adjacent to 
a more evolved young stellar object. Similarly, in ISOSS 
J20153+3453 one massive clump of about 100 M is present 
and a cluster of sources is detected in the vicinity containing 
a presumably deeply embedded accreting protostar with a 
current mass of ~ 1.5 M Q . Its evolutionary state corresponds 
to Class 0/1 and there is evidence for outflow activity from 
this object. 

6. The source SMM1 in ISOSS J20298+3559 probably repre- 
sents a young accreting intermediate-mass star precursor em- 
bedded in a molecular clump of ~ 1OM . Towards SMM4 
in the same region we identify an embedded proto- or pre- 
main-sequence star of 3.2 - 4.9 M that likely evolved from 
the associated clump of ~ 60 M . 

7. In the region ISOSS J22478+6357 a candidate protostar of 
6-8.5 M is found that is embedded in the associated SMM1 
East clump (~ 60 M ). 

8. We do not detect stellar precursors that could have cur- 
rent masses of 1OM or more. Therefore we cannot be 
sure whether such objects will emerge within the identified 
clumps. However, the presence of intermediate-mass proto- 
or pre-main-sequence star candidates and the large clump 
masses indicate that one can expect that also high-mass stars 
are formed in these systems. 

9. Our study reveals that the search for regions containing large 
fractions of cold material in the far-infrared has success- 
fully identified star-forming regions associated with cold and 
massive clumps. Compared to sources selected at shorter 
wavelengths they may represent less evolved stages of high- 
mass star-formation and in that respect substantial impact 
can be expected from future (space) missions like Herschel 
and Planck. Furthermore this study demonstrates that the 
detailed star-forming content can only be accessed with a 
multi-wavelength approach including sensitive mid- and far- 
infrared observations. 
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